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Abstract. The equations which determine the response of a charged particle mov- 

I ing in a magnetic field to an incident gravitational wave(GW) are derived in the 

' linearized approximation to general relativity. We briefly discuss several astro- 

, physical applications of the derived formulae taking into account the resonance 

, between the wave and the particle's motion which occurs at uig = 2f2, whenever 

I , the GW is parallel to the constant magnetic field. In the case where the GW 

I^J^, is perpendicular to the constant magnetic field, magnetic resonances appear at 

^ , LUg — fl and u)g — 2Q. Such resonant mechanism may be useful to build models of 

•1— I . 

. GW driven cyclotron emitters. 
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1. Introduction 

Papadopoulos and Esposito (1981) discussed the perturbations of the Larmor orbits 
in the presence of a gravitational wave(GW) and estimated the consequent magnetic 
bremsstrahlung. They has shown that it is possible to identify the presence of GW in 
macroscopic systems by detecting the shifts in the spectrum of the electromagnetic ra- 
diation (cyclotron) given off by charged particles as they interact with a GW. 

Recently, the motion of a relativistic charged particle in a constant magnetic field per- 
turbed by GW incident along the direction of the magnetic field has been examined(Jan- 
Willem van Holten 1999 and references therein). In the same work, a generalized energy 
conservation law to compute the variations of the kinetic energy of the particle during 
the passage of the GW has been derived and explicit computations in the orbit of the 
charged particle due to the GW has been obtained. 
Send offprint requests to: D. Papadopoulos 
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In this paper, we discuss the interaction of GW with a charged gyrating particle in the 
presence of a constant magnetic field across the z-axis in the frame of linearized theory 
of gravity. 

When the GW propagates parallel to the magnetic field with a frequency iVg, the 
coupling of a gyrating particle with the GW becomes very strong at the resonance which 
occurs between the GW and the Larmor orbits. The resonance is at twice the Larmor 
frequency(n = ^) e.g ujg = 20 (Macedo and Nelson (1990)). 

In the case that the GW propagates perpendicular to the magnetic field, the interac- 
tion again becomes extremely efficient at the resonances, w = Q and co = 2U. 

In both cases we verify that close to gyro resonances, the obtained spectrum of the 
produced cyclotron radiation, becomes comparable to the spectrum of the initially gy- 
rating particle especially in the vicinity of a source producing the GW . 

Our results suggest that a) the linear theory breaks down at the resonances and the 
interaction of the GW with gyrating particles becomes very efficient and, b) even in the 
linear theory, support the discrepancy on the estimations for the cyclotron damping radi- 
ation recently discussed by M. Servin (Martin Servin et al (2001)) and Kleidis (K. Kleidis 
et al (1996); K. Kleidis ct al (1995)), since in Kleidis work, the problem is examined in 
the non-linear theory where magnetic resonaces occur and some of them are overlapped. 

The paper organized as following. In sec. II, we derive the equations of motion in 
the linearized theory. In sec. II wc discuss the interaction of the gravitational wave to 
the magnetic field when the GW is parallel to the magnetic field. In sec. IV we discuss 
the same problem assuming that the GW is perpendicular to the magnetic field. The 
obtained results are discussed in sec.V 
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2. Derivation of the equations of motion 

In the linearized approximation to general relativity the metric tensor is decomposed in 
the fashion 

9i3 = Vij + hij (1) 

where the elements hij are small compared to unity. By imposing the condition 



iH-^if^^j=o (2) 

we reduce the vacuum field equations to homogeneous wave equations for all com- 
ponents of . The gravitational field is then described by a symmetric traceless, diver- 
genceless tensor with two independent space components. Thus, the square of the line 
element is 

ds^ = {rjij + hij)dxUx^ = {dx°f - {dx'^f + Kpdx'^dx^ (3) 
where Greek indices take values 1,2,3 and Latin 0,1,2,3. 

The components of the covariant four- velocity, consistent with the linearized theory, 

are 

dx'^ 1 



where are the components of the four- velocity and the same quantities with 

the subscript M distinguish the special-relativistic Minkowski values. 

The equations of motion of a test particle, which is moving in the presence of an 
electromagnetic field in the space-time defined by Eq.(l) are given by 



^+T%uiu'^ = ^F^^uu (6) 
ds ■' mc^ 

where the right-hand side is the inhomogeneous driving term determined by the elec- 
tromagnetic field in the space- time defined by the Eq.(l) 
For the metric (1), the nonzero Christoffel symbols are 



■pO "^1. pa La 



^, = ^{hl^ + hl0-h-p^) (7) 



From Eqs.(6),(4),(5) and (7), the equations of motions take the form: 



4 
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and 
where 



Finally, the equations of motion (8) and (9), in the Newtonian and linearized limit, 
reduce to the equations: 

^ + v''hjk,tv^ + lv''[hji,k + hki,j - hjk,i]v^v>' = ^(r?'« - (11) 

To make further progress with the derived equations of motion (11), we consider the 
gravitational wave which is characterized by the wave vector 

fc" = -(sin(6'),0,cos(6')), and = Wka (12) 
and one of two possible states of polarization given by 



hij = ho{e\e] - eh]) exp [i'^ix^ sm(e) + cos(e) - ct)] (13) 

where ho is the amplitude of the gravitational wave and Ug = 2'KUg is the angular 
frequency of the GW. 

The vectors and have space components only and satisfy the conditions 



gi'^ej, = e^'^e^ = 1, with, k^'el = k^'el = (14) 
Conditions (14) imply 

=(cos(^),0,-sin(^)), =(0,1,0) (15) 

Under the above consideration wo proceed to the following two cases, a) the GW is 
parallel to the magnetic field and b) the GW is perpendicular to the magnetic field which 
will be analyzed in the following two section. 
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3. The GW is parallel to the magnetic field 

We choose the electromagnetic field to be 



5 



F- ■ 
^ 1,1 





-iJs 

ils 





(16) 



where H'^ = H = constant. 

Wc choose the gravitational wave to propagate parallel to the magnetic field e.g., in 
Eq.(13) we obtain 6 — Q and thus h = hu = — /122 = hoexp{^{z — ct)). Eq.(ll) with 
the aid of Eq.(16) yields: 



dv 



dt ^ ^ c ^ 



(17) 



— + f2t; = hiu}gV^[l 

dt " ^ c 



(18) 



1 



(19) 



dt "■"2c-^^^ 

where Vt= To solve the system of Eqs. (17-19), we decompose the components of 
the 3-velocity as follows: 



(20) 



where the subscript zero means zero order in the sense that ho = 0, while the subscript 
one means first order in the sense that ho ^0. 

The perturbed equations of motion are derived from Eqs. (17-19) and (20). Thus, 
after some straightforward calculations (see Appendix B) we find the solution(Macedo 
and Nelson (1990)): 



0:^ VoTCOs{flt + a) + hoVoTj^ — ^^^^{cos[kgZ + {'^ — ujg)t] — COS {kgZ — ^t)} (21) 

(2i2 — ^g) 



-voT sm{ilt + a) + hoVoT 



-^^^^{sin [kgz + u;g)t] - sin {kgZ - Qt)} (22) 



^3 ^ '^^^),^^if,^,){ ~. 



ho Vq 
2 c 

+ LOg) 



exp (i(2f2 — LOg)t 



exp {02^1 + ujg)t\} 



(23) 
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where Vqj: = ''Jx ~^ "^y ~ constant and a = constant. 

It is evident that in Eqs. (21-23) if ho = 0, we obtain the components of the space 
velocity of the initial gyrating charged particle. If ho / 0, the Eqs. (21-23) reveal that 
the gyrating charged particle diverts from its initial plane orbit and moves into a helical 
trajectory. Now the vector v = {v^,v'^,v^) does not move in a circle, but on the surface 
of cone with its axis along the H. From Eqs. (21-23) we conclude that the particle is 
accelerated at the resonace ojg — 217. The existence of the resonace at uig = 2f2, in 
Eqs. (21-23) is due to GW. Because of the resonance the gyrating particle gains kinetic 
energy. Thus, if Eq and Ei are the kinetic energy of the particle before and after the 
interaction with the GW, respectively; the energy gained by the particle in a period, let 
say T, is given by the average of the ration ^ e.g. 



1=1- r ^dt ^l- 2ho^-^ cos {kgz) (24) 
Obviously, as ojg approaches 20 taking value between fl and 2f2, the factor, ^""^^ 



becomes negative making the term multiplied by ho, to approach to plus infinity with 
positive values. This suggests extra emission of cyclotron radiation which will change the 
spectra distribution of the radiation of the initial gyrating charged particle and transfer 
of energy from the GW to the particle. 

Integrating the Eqs. (21-23) we obtain the parametric equations of motion of the 
charged gyrating particle interacting with a GW in the presence of a constant magnetic 
field across the z-axis. These are 



^ m _ ^ , ^ ^-^g , sm[kgz + {n-ujg)t] sin {kgZ-ftt) ^ 
^(i)(i) - + f'oVoT^^^i ^^3- ^ } (25) 



Q-tOg COS [kg z + {n - (jjg)t] COS {kgz - nt) 
y(i)(i) = 2/(01) - hoVoT^^^i + ^ } (26) 

and 

Z(i){t) = ^(01) + Y^0T{—)exp{ikgz){ ^^2_^2 

+ l^W^ " ^ WT^ ^ "^^'^^ ^''^ 

where a:(oi) ,2/(oi) and Z(oi) are constants of integration. 

The intensity of the radiation per solid angle per unit interval of frequency produced 
from a charge test particle moving in the presence of a magnetic field which interacts 
with the GW, maybe obtained from the relation(Landau 1975): 

- 4Srl £*»p(i.'(*-=^))[„x(„xb)ir (28) 

where lo' is the frequency of the outgoing radiation, n = sin {6)i + cos ((?)k, R is a vector 
which joins the charged particle with the observer, b is the velocity of the charge particle. 
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We carry out the integral of Eq.(28) neglecting terms of the order Thus, we find: 

,o ^ 9 o oo n 



— oo 

- 4/10 cot^ W cos Jf ($)} (29) 

where ro = ;n[l + /i0 2S7^S^^ sin (A;g0)] =;f7[l + 0(l/c2)], $ = sin (0), J, is the 

Bessel function of first kind and J; its first ordinary derivative in terms of its argument. 
For simplicity we call 

L = [cot' {9) JH^) + ^{j'irm (30) 
and 

T = -4ho cot' (6) cos {kgz) 2r'^' Jfi^) (31) 

It is evident that, as cjg approaches 2f7+, the factor 2n-J ^ °° making the term 
T tend to minus infinity. But, If ujg approaches 2Q taking value between Q and 2il, the 
factor, ^2^u! ' becomes negative making the term T, to approach to plus infinity with 
positive values. Also, from Eq.(31), wc sec that the divergence of the term T is faster 
as we approach to the source producing GW. Nevertheless, in the linearized theory of 
gravitation, we have to approach the resonant in such a way that the term T remains 
always below the term L, otherwise the linear theory breaks down. 



4. The GW is perpendicular to the magnetic field 

We shall now consider the case where the GW propagates perpendicular to the unper- 
turbed magnetic field = (0, 0, H^) = const. Thus, in a reference frame where has 
the z-direction (as for the previous case), but the GW propagates along the x-direction, 
the two non-vanishing components of the GW are given by Eqs.(13) setting 6 = 7r/2 and 

/i33 = —h22 — ho exp — (x - ct) (32) 
c 

Subsequently, from Eqs.(ll),(20) and (32) we obtain the following equations of motion: 
= Vq+vI= vqt cos {fit + a) + voTho{C cos {fit) 

— [Acos{kgZ — ojgt) — Bsm{kgZ — ujgt)]} (33) 



-uqt sin {fit + a) — 'L'ot^o{C' sin {fit) 



+ [Asni{kgZ — UJgt) — B C0S{kgZ — UJgt)]} 



(34) 



and 



= vl + vl = a-a{l )/l22 



(35) 
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where, a = constant and may be chosen equal zero, 



C = — ^1— [i^^^) ay 9J . (3g^ 



COS (lit + 2a) 



fl — UJg — UJg 

voT _ ^.(3» + ^.) (20i + 2a)] (37) 

and 

B = sin (nt + 2a) - ^ — r sin (20t + 2a) (38) 

Obviously, if in Eqs. (33-35) ho = 0, then we obtain the space velocities of the initial 
gyrating charged particle. 

If ho ^ 0, then we have gyrating motion again, but magnetic resonances appear at 
ojg = ±Q,,u)g = 20 and ojg = 3Cl. 

The existence of the above resonances, in Eqs. (33-35) is due to GW. Because of those 
resonances the gyrating particle gains kinetic energy. Thus, as in section (III) we verify 
that the energy gained by the particle in a period, let say T, is given averaging the ration 
l^e.g. 



()- 



^2 = ;^ / ^dt f« 1 + 2hoC « 1 + 2hoj^ '—- r (39) 

J- Jo ^0 



— rU ^. 1 -I- ')h^r a-. 1 -I- 7hr, 

{n - (Jg){2n - UJg) 

Obviously, if Wg approaches 20 from the right or O from the left, the factor, 
(Q-u}g)(2n-ojg) ' becomes positive making the term multiplied by ho, to approach to plus 
infinity with positive values suggesting changes to the spectra distribution of the radiation 
of the initial gyrating charged particle. 

Integrating Eqs. (33-35) we derive the parametric equations of motion which are: 

a:{t) = ^ sin (nt) + hoVorC-^^^C + X^] (40) 
c il 



y^t) = ^ COS (Ot) + hoVoA^-^^C + Y,] (41) 

C O 



z{t) = ut + i—{l--)h22 (42) 

Ulg C 

where, a = constant and may be chosen equal zero, the expression for and are 
given explicitly in the Appendix A. 

In the Eqs. (40) and (41) a drift term, with resonances at Ug = fl^ujg ~ 20, is present. 
This drift term can generate electric currents. Those currents arc sources of secondary 
electromagnetic waves. For further details see (Macedo and Nelson (1990)). 
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Following the same procedure as in in section III, Eq.(28) reduce 



d^ladoj' c ^ 6 



-oo 

■)2 



.0 COS (2.) cot^ (43) 
As in sec. Ill, we call 

L = [coeieym + '^iJ'ifm (44) 

and 

= -/.ocos(2.)cot-W(^) ^^_j;;;_^^^ jf(^) (45) 

Now the term T2 has two magnetic resonances. It is evident that, as Wg approaches 
2f]~, or 0+ and 6 G (0, f ), the term T2 tend to plus infinity. Also, from Eq.(45), we see 
that the divergence of the term T2 becomes faster as we approach to the source producing 
GW and as we are dealing with ultra relativistic particles where ratio takes higher 
values. Again, as in in the section (III), in the linearized theory of gravitation, we have 
to approach to the resonances carefully!, in the sense that the term T2 should not exceed 
term L, otherwise the theory breaks down. 
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5. Conclusions 

In this article we pose the foUowing problem. If /iq = 0, we reproduce a well known formula 
for the spectrum distribution of a gyrating charged particle in both cases e.g when the GW 
is parallel and perpendicular to the magnetic field. In this case the angular distribution 
of the gyro radiation is highly anisotropic. The radiation is concentrated mainly in the 
plane of the orbit. 

If /lo 7^ 0, we are dealing with the interaction of a GW with a gyrating charged 
particle. We have distinguished two cases: 

(a) . The GW propagates parallel to the constant magnetic field: 

Because of the GW, the gyrating charged particle diverts from its initial plane orbit 
and starts to move across a helical trajectory. Also, because of the GW, the spectnun 
of the produced cyclotron radiation isolated by a factor proportional to ho, namely T, 
in which a resonant at ujg ~ 2fl appears. If ujg approaches 2il taking values in the 
interval Ii = (fi, 217), the term T, approaches to plus infinity indicating that the existence 
of the resonant interaction between the charged particle and the GW, can lead to a 
strong emission of cyclotron radiation, even in the linear theory of gravity. The suggested 
mechanism of cyclotron radiation could be useful to the astrophysicists especially in the 
case they make simultaneously observations of the so obtained cyclotron radiation and the 
GW. Therefore, knowing that astrophysicists are looking to detect GW at frequencies 
between i/g = (10^ — lQp)Hz (Cutler, Thorne 2002), magnetic resonance may occur 
whenever the magnetic fields could be between H k, {1. — 2.)G [vg = IQ'^Hz), H f« 
(10 — 21)G (ug = IQ^Hz) and in both cases the T term becomes positive and comparable 
to the L term . But, if we want to use electromagnetic radiation for indirect detection of 
GW, then other frequencies are also important, e.g. lower frequencies of Ug = {1 — 10)Hz 
coming from binary neutron stars before coalescence, the corresponding magnetic fields 
are H « (0.01 - 0.02)0 {vg = IHz), H « (0.1 - 0.2)G {vg = 10Hz)or high frequencies 
fg = (10^ — 10'^)Hz, due to supernovae explosions, to normal modes of pulsating neutron 
stars or stellar size black holes. In the latter case the corresponding magnetic fields are 
« (107-214)G. 

Outside interval Ii the particle seems to lose energy due to destructive interference 
of the two oscillators e.g. the sinusoidal gravitational wave and the gyrating particle. In 
this case, the actual loss described by the negative term T is of the order Hq. 

(b) . The GW propagates in the x-direction e.g perpendicular to the constant magnetic 
field. 

In this case, we verify that, 1) the gyrating particle remains on the initial plain of 
orbit described by Eqs. (33-35) and for a certain 6 = Oq, let's say, £ (0, j), at ujg — 
and LUg = 2Q the corresponding term T2 diverges. 2) If ujg approaches to Q from the 
left, or to 2Q from the right, whereas G (f i f )> the term T2 becomes positive, and 
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approaches to plus infinity, unless 9o G (0, ^). For the frequencies mentioned above, two 
magnetic resonaces may occur for the same values of ujg at the same range of the magnetic 
field. But if tOg approaches to one or the other resonance taking values in the interval /i, 
whereas G (f j f )j the term T2 becomes negative (unless G (0, f )), indicating that 
the particle lose energy due to the same reason mentioned in paragraph (a). Another 
interesting feature is that the term T2 is proportional to the ratio , which means that 
high relativist particles support the term T2 to become more significant. Nevertheless, in 
both cases, the strength of the cyclotron radiation described by the ratios and 
remain constant for large values of l(large frequencies). This may seen in figure 1 where 
we plot the above mentioned ratios for 6q — 80°, obtaining the amplitude of the GW to 
be ho — 10"^^ and its frequency — lO^^. 



l.OE-5 



T/Q 



^ l.OE-10 - 



o 

t - 

I ' 

& : 



i.oE-12 ^ — , — ^ — , — ^ — , — ^ — , — ^ — , — 

20 40 60 80 100 

1-modes 

Fig. 1. The strengths of cyclotron radiation in the cases where H is parallel to the 
A;(T/L1) and H is perpendicular to the k 



However, we have to point out that, in both cases, the terms multiplied by Hq should 
not exceed the term L, otherwise the linearized theory of gravity breaks down. 

In the nonlinear theory, the problem could be more interesting. The problem some- 
how has been examined from a dynamical point of view, considering the equations of 
motion from a Hamiltonian, ( Varvoglis and Papadopoulos (1992), Kleidis Varvoglis and 
Papadopoulos (1993), Kleidis Varvoglis, Papadopoulos and Esposito (1995) and Kleidis 
Varvoglis and Papadopoulos(1996)), and integrating them numerically. In this case, the 
interaction of the GW with gyrating charged particle exhibits resonances and in several 
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cases chaotic behavior. We intend to discuss the problem in the non Unear theory in a 
forthcoming paper. 
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Appendix A 

In the Eqs. (33-35) theXh and are: 



13 



sin [fc.^ + - + [^y^ sin - (O + a;,)i] 



and 



+ TT-^^ — sin [fc„a; + (O - w„;„jj , 

Cl?^ 1 1 

+ 7T77^^ tt^tf: T sin \k„x - {20. + ijj„)t] + — i 

n n{Q.-2ujg) 2nujl 



- sin [fcgX + (20 — Wg)^]] 



4c ^n-Wg (20-|-a;g)(3fi 



~y/\-- ~gy . -g/v - 

a;^(3n + u;g) I 

(3fi - a;g)(4n2 _ u^f) (fi + a;g)(3fi 



= — 1— ttH 7;^ r cos (Kga; — Wof) — -— — r- 

^ O a;g(0-a;g) ^' 2(20 -Wg)^ 

O — O^g COS (fcgX — (f2 — Wg)*) COS ( /Cg " ( + Wg ) t ) . 

2(20 -Wg)^ O-Wg O + Wg 



O .COS (/CgX — (O — Wg)^) COS (/CgX — (O + 



O — Wg 



2(20 - Wg) ^ O - Wg 

— COS iknx) \ — - —— 

^ ^ ^^a;g(0-Wg) (20 



-wg) (02-^2) (20-u;g)(0+a;g)^ 

— ojgt) sin (20i) 



2w, 



c 

T~l7?; N *^°S {kgX - UJgt) + — : 

4c (0-Wg) (0 - tJg 

Wg(30 + Wg) COS (fcgX — (20 — Wg)^) 

^1 "0_,., + 20 + UJg 



Wg(30 + Wg) COS (fcgX — (20 — Wg)^) COS ( ^ 

2(0 + Wg)(30-a;g)^ 20 - Wg ' 

...^ 20a;g(30 + a;g) 

cos [l^gXn^ ^ ^^^^3^ _ ^^)(4Q2 _ ^2)J> 



-sin{knX — UJnL 

,)t) cos {kgX — (20 + iOg)t) 



g)(30- 
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Appendix B 

We start with the Eqs.(17-19) 



^-Qv^ = -hiLOgV^[l - —] (48) 



— + = hiLOgV^[l - —] (49) 



^ = -^hj-^^,[{v'r - iv^] (50) 
We write the Eqs.(48),(49) as follows: 

^^'''^''''^ + n[v^ + iv^] = hiiogiv^ - iv^)[l - ^] (51) 

In order to solve the Eq.(51), we decompose the components of the 3- velocity as 
follows: 

v^^vl+v\, v'^^vl + vl, w^-O + i;? (52) 

where the subscript zero means zero order in the sense that /iq = 0, while the subscript 
one means first order in the sense that /iq ^ 0. 

We substitute Eqs.(52) into Eq.(51) and the perturbed equation now reads : 



d[{vl + v\) + ijvl + ivl)] 

at 



+ in[{vl + vl) + i{vl + vl'' 



= hiujg[{vl + vl) - i{vl + vl][l - -j] (53) 



I. Zero Order Equations 



This yields 



■^0 = ^OT cos{ilt + a) and Vq = —vqt sm{Qt + a) (55) 
where a = constan and Vqj, = v^+Vy = constan. 
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II. First Order equations 
From Eq.(53), we have 

^^''^'^''''^ + m{vl + ivf) = hiojgivl - ivl) (56) 
Notice that on right hand side, the factor h{vQ — iVq) gives 

h{vl - ivl) = /lowore'^'^''''"'^''*^ [cos {Qt + a)-i sin {Qt + a)] 

= /ioUoTe'('=''^-'"''*)e'("*+") (57) 

Now from Eqs. (56) and (57) we have 



We treat Eq.(58) as an ordinary first order differential equation with the initial con- 
ditions, if f = then {vl + ivf) = 0. Thus, we have 



g 



Upon the consideration of the initial conditions we have 
C = -/ioi;oT7T7^^^e'('=«^+«) (60) 
Eventually, from the Eqs. (59) and (60) we obtain 



(y\ + ivl) = hoVoT ^Jr" e'(''''^+°)[e^*("-"'') - e"**"] (61) 
or 



1 '^Q 

Vl = hovoT — {cos [kg z + {fl — LOg)t] — COS {kg z — Qt)} (62) 

2iZ UJg 



vl = hoVoT ^f^^ — {sm[kgZ + - ujg)t\- s\n.{kgZ -Vlt)] (63) 

Z\l LOg 

Furthermore, following the same method, from Eq.(19) we find 



2 



V 



3 _ ho ff'm\ w 



1 = — uot( — )oxp(ifcgz){ 



2 ""^^ C ^""^^"^3"^^4f]2_^2 

[ (20^ exp i^i2n - .g)t - exp (^(2^^ + u;g)t]} (64) 
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